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Abstract 

The GacS/GacA two-component regulatory system activates the production of 
secondary metabolites including phenazines crucial for biological control activ- 
ity in Pseudomonas chlororaphis 30-84. To better understand the role of the Gac 
system on phenazine regulation, transcriptomic analyses were conducted by 
comparing the wild-type strain to a gacA mutant. RNA-seq analysis identified 
771 genes under GacA control, including many novel genes. Consistent with 
previous findings, phenazine biosynthetic genes were significantly downregulat- 
ed in a gacA mutant. The transcript abundances of phenazine regulatory genes 
such as phzl, phzR, iopA, iopB, rpoS, and pip also were reduced. Moreover, the 
transcript abundance of three noncoding RNAs (ncRNAs) including rsmX, 
rsmY, and rsmZ was significantly decreased by gacA mutation consistent with 
the presence of consensus GacA-binding sites associated with their promoters. 
Our results also demonstrated that constitutive expression of rsmZ from a non- 
gac regulated promoter resulted in complete restoration of AT-acyl-homoserine 
lactone (AHL) and phenazine production as well as the expression of other gac- 
dependent secondary metabolites in gac mutants. The role of RsmA and RsmE 
in phenazine production also was investigated. Overexpression of rsmE, but not 
rsmA, resulted in decreased AHL and phenazine production in P. chlororaphis, 
and only a mutation in rsmE bypassed the requirement for GacA in phenazine 
gene expression. In contrast, constitutive expression of the phzI/phzR quorum 
sensing system did not rescue phenazine production in the gacA mutant, indi- 
cating the direct posttranscriptional control by Gac on the phenazine biosyn- 
thetic genes. On the basis of these results, we propose a model to illustrate the 
hierarchic role of phenazine regulators modulated by Gac in the control of 
phenazine production. The transcriptomic analysis also was used to identify 
additional genes regulated by GacA that may contribute to the biological 
control capability of strain 30-84. 



Introduction 

The GacS/GacA two-component signal transduction sys- 
tem (TCST) is highly conserved among Gram-negative 
bacteria, including many beneficial biological control and 
plant pathogenic bacteria (Heeb and Haas 2001). The Gac 
system was originally identified as a global activator of 



antibiotic and cyanide synthesis in Pseudomonas, hence 
the acronym Gac (Haas and Keel 2003). It is now known 
that Gac regulates the expression of secondary metabolites 
(including hydrogen cyanide [HCN], phenazines, 2,4- 
diacetylphloroglucinol [DAPG], pyoluteorin, and the phy- 
tohormone indole-3-acetic acid), extracellular enzymes, 
and several carbon storage compounds as well as an oxi- 
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dative stress response and other functions that play roles 
in biological control or plant pathogenicity (Hassan et al. 
2010). 

Global regulation by Gac in many Pseudomonas species 
was shown to act through the small RNA-binding 
proteins RsmA and RsmE, members of the RsmA/CsrA 
family (Reimmann et al. 2005; Lapouge et al. 2008). 
According to the current model, these small RNA-binding 
proteins act as negative regulators of gene expression by 
binding to 5'-GGA-3' motifs located in the 5' leading 
sequence of target mRNAs (Blumer et al. 1999; Heeb and 
Haas 2001). The binding of these proteins reduces target 
protein levels either by a reduction in mRNA translation 
(by blocking the recruitment of the 30S ribosomal sub- 
unit), by a reduction in mRNA stability (by targeting 
messages for degradation), or both (Dubey et al. 2003). 
Gac also is required for the expression of the small regu- 
latory noncoding RNAs (ncRNAs), rsmX, rsmY, and rsmZ 
(Lapouge et al. 2008; Humair et al. 2010). These ncRNAs 
have multiple stemloop structures that interact with 
RsmA and RsmE, sequestering them, and relieving post- 
transcriptional repression (Babitzke and Romeo 2007). 

The rhizosphere-colonizing bacterium Pseudomonas 
chlororaphis 30-84 is a member of a group of phenazine- 
producing bacteria with beneficial agronomic applications 
(Pierson and Pierson 2010). Strain 30-84 is a biological 
control agent capable of suppressing take-all disease of 
wheat caused by the fungal pathogen Gaeumannomyces 
graminis var. tritici (Ggt) (Pierson and Thomashow 1992). 
Strain 30-84 produces phenazine-l-carboxylic acid (PCA), 
2-hydroxy-phenazine-l-carboxylic acid (20HPCA), and a 
small amount of 2-hydroxy-phenazine (20HPZ) (Pierson 
and Thomashow 1992). Phenazines are responsible for 
the majority of strain 30-84's ability to inhibit fungal 
pathogens and have been shown to be important for 
30-84 to persist in the rhizosphere (Mazzola et al. 1992) 
and to form biofilms in vitro (Maddula et al. 2006). 

Phenazine biosynthesis in 30-84 is controlled by a 
complex regulatory network involving quorum sensing, 
Gac and other TCSTs, sigma factor rpoS, and other regula- 
tory genes (Pierson and Pierson 2010). Phenazine biosyn- 
thesis is regulated directly by the PhzR/PhzI quorum 
sensing system (Pierson et al. 1994; Wood and Pierson 
1996; Wood et al. 1997). PhzR is a LuxR family transcrip- 
tional regulator that activates the expression of the phen- 
azine biosynthetic genes in response to the accumulation of 
N-acyl-homoserine lactone (AHL) signals. Mutation in 
phzR abolishes the expression of the phenazine biosynthetic 
genes and multiple copies of phzR in trans significantly 
increases phenazine production (Pierson et al. 1994). Phzl, 
a Luxl homolog, is an AHL synthase that produces the AHL 
signals to which PhzR is most responsive (Wood et al. 
1997). In 30-84, mutations in Gac result in loss of phen- 



azine, HCN, exoprotease, lipase, gelatinase, and AHL pro- 
duction as well as losses in the capacity for pathogen 
inhibition, biofilm formation, and rhizosphere competence 
(Chancey et al. 1999; Maddula et al. 2006). Other genes 
with demonstrated roles in phenazine production in P. 
chlororaphis include pip (pjienazine inducing protein), the 
rpeA/rpeB TCST, and iopA/iopB (inducers of phenazine) 
(Girard et al. 2006a; van Rij 2006; Wang et al. 2012b). Pip, 
a homolog of the TetR family of transcriptional regulators 
(Ramos et al. 2005) was shown to promote phenazine pro- 
duction in 30-84 and P. chlororaphis PCL1391 by enhancing 
phzl and phzR expression (Girard et al. 2006a; Wang et al. 
2012b). In both P. chlororaphis 30-84 and PCL1391, the 
expression of pip is regulated by the sigma factor rpoS. In 
30-84, the RpeA/RpeB TCST also regulates pip (Wang et al. 
2012b). In PCL1391, two additional genes, iopA and iopB, 
positively contribute to phenazine production through 
RpoS (van Rij 2006). Pseudomonas chlororaphis 30-84 also 
has a second quorum sensing system, Csal/CsaR, that con- 
trols exoprotease production and cell aggregation, but plays 
a minor role in phenazine regulation under laboratory con- 
ditions (Zhang and Pierson 2001). 

GacS/GacA also negatively regulates a spectrum of traits 
by mechanisms that are as yet largely uncharacterized. For 
example, inactivation of GacA resulted in the differential 
expression of 635/6,147 (10%) of the genes in P. protegens 
(fluorescens) strain Pf-5 (Hassan et al. 2010). Of these, 288 
genes were upregulated in a gacA mutant. On solid med- 
ium colonies of gac mutants often appear hyper-fluores- 
cent and larger than wild type (WT), suggesting that traits 
related to iron acquisition (e.g., siderophores) or motility 
may be among the traits negatively regulated by Gac 
(Chancey et al. 1999). Negative regulation by Gac is par- 
ticularly interesting given that many Pseudomonas species 
identified for their biological control potential exhibit 
phenotypic variation resulting from spontaneous mutation 
in gacS or gacA (van den Broek et al. 2005). Spontaneous 
gac mutants are especially prominent in fermentation cul- 
ture, where they often outgrow the WT population result- 
ing in a deficiency of secondary metabolites essential for 
biological control activity (Duffy and Defago 2000). In 
comparison, the prevalence of gac mutants in natural envi- 
ronments (e.g., in the rhizosphere) suggests some benefit 
to the Gac~ phenotype and thus selection at some level to 
maintain it. In contrast to liquid culture, in the rhizo- 
sphere both the WT (Gac + phenotype) and the Gac 
mutants (Gac - phenotype) survive better in mixed popu- 
lations than in uniform populations, demonstrating a ben- 
efit to maintaining phenotypic variation (Chancey et al. 
2002). More recently, it was shown that the WT and gac 
mutants interact mutualistically in biofilms (Driscoll et al., 
2011). The improved survival of both Gac + and Gac~ phe- 
notypes in mixed populations suggests that mutants 
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arising de novo or introduced within the inoculum 
contribute to the survival of 30-84 in the rhizosphere, 
however, the molecular mechanisms responsible remain 
unknown. 

In this study, the functionality of major components of 
the Gac/Rsm signal transduction pathway in 30-84 was 
demonstrated. Transcriptomic analysis using RNA-seq 
was performed comparing WT and a gacA mutant to 
examine whether the expression of known phenazine reg- 
ulatory genes including phzR/phzI, pip, rpoS, and rpeAl 
rpeB were controlled by GacA. We used mutation and 
complementation analysis of genes within the Gac/Rsm 
pathway as well as phenazine regulatory genes to deter- 
mine their hierarchic role in the control of phenazine 
production. This study expands the current knowledge of 
the Gac/Rsm regulon, especially its role in phenazine reg- 
ulation, and based on our results a model describing the 
regulatory network controlling phenazine gene expression 
in 30-84 is proposed. The transcriptomic analysis also 
facilitated the identification of additional genes regulated 
by GacA that may contribute to the biological control 
capability of 30-84. Of special interest were genes with 
increased expression in the gacA mutant as these may 
provide insights into the benefits of the Gac~ phenotype 
observed in mixed populations. 

Experimental Procedures 

Bacterial stains and growth conditions 

Bacterial strains and plasmids are listed in Table SI. Oli- 
gonucleotides and primers are listed in Table S2. Liquid 
Luria Broth (LB) medium or AB minimal medium sup- 
plemented with 2% casamino acids (AB + 2% CAA) (Dif- 
co, Becton Dickinson and Company, Franklin Lakes, NJ) 
were used for culturing P. chlororaphis as described previ- 
ously (Wang et al. 2012b). The following antibiotics were 
added to the medium when necessary: ampicillin (Ap) 
100 /(g/mL, kanamycin (Km) 50 /(g/mL, rifampin (Rif) 
100 /(g/mL, tetracycline (Tc) 50 /(g/mL, and gentamicin 
30 /(g/mL. 

DNA manipulation and sequence analysis 

Standard procedures were used for plasmid isolation, 
cloning, restriction enzyme digestion, and T4 DNA 
ligation (Sambrook and Russell 2001). Polymerase chain 
reaction (PCR) was carried out using Invitrogen Taq 
DNA polymerase (Life Technologies, Carlsbad, CA) at 
95°C for 5 min, followed by 30 cycles of 95°C for 30 sec, 
60°C for 30 sec, and 72°C for 90 sec, and a final elonga- 
tion step of 70°C for 10 min. DNA sequencing was 
performed at the Genome Technology Lab within the 
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Texas A&M University Institute for Plant Genomics and 
Biotechnology. Nucleotide and amino acid homology 
searches were conducted using the blast programs 
(BLASTN and BLASTP, respectively) on the National 
Center for Biotechnology Information (NCBI) website 
(http://www.ncbi.nlm.nih.gov/BLAST). The P. chlororaphis 
subspecies aureofaciens 30-84 genome sequence was 
recently deposited in GenBank (GenBank Accession #: 
PRJNA67533). The P. chlororaphis 30-84 phenazine bio- 
synthetic operon contains seven conserved biosynthetic 
genes phzXYFABCD (Pierson et al. 1995; Mavrodi et al. 
1998), which correspond to phzABCDEFG (according to 
the P. fluorescens nomenclature, Mavrodi et al. 1998). 
Here, we use the P. chlororaphis nomenclature to conform 
to the original literature. 

RNA preparation for RNA-seq and qPCR 
analyses 

Three biological replicates of every strain were started from 
single colonies located on three separate plates containing 
AB + 2% CAA and then transferred to 10 mL AB + 2% 
CAA broth. All cultures were grown at 37°C with shaking 
(150 rpm) to an approximate OD M0 = 1-2. Cell cultures 
collected at OD 600 = 1.2 were diluted to OD 600 = 0.3 with 
AB + 2% CAA broth. RNA extraction was performed as 
described previously (Wang et al. 2012b) with one excep- 
tion: contaminating genomic DNA was removed off-col- 
umn with Turbo DNA-free DNAse (Life Technologies, 
Carlsbad, CA). Elimination of contaminating DNA was 
confirmed via qPCR amplification of the rpoD gene with 
SYBR green® dye on an ABI 9400HT PCR machine (Life 
Technologies, Carlsbad, CA). RNA samples were ethanol- 
precipitated and resuspended in 0.1% diethylpyrocarbonate 
(DEPC). Ribosomal RNA (rRNA) was depleted from 
~9 \i% of total RNA using the RiboZero rRNA depletion kit 
(for Gram-negative bacteria, Epicentre Biotechnologies, 
Madison, WI). Two separate samples were prepared for 
each treatment using this protocol. At the elution step, the 
two samples were pooled and concentrated using RiboMi- 
nus concentration modules (Life Technologies, Carlsbad, 
CA). RNA quantification throughout was achieved using a 
GE NanoVue Plus spectrophotometer (GE Healthcare Bio- 
Sciences Corp., Piscataway, NJ) and RNA quality was mon- 
itored during preparation for library construction and dur- 
ing library construction with an Agilent 2100 Bioanalyzer 
(Agilent Technologies, Santa Clara, CA) at the Texas A&M 
GTL. 

Library construction and sequencing 

Strand-specific cDNA libraries were constructed using the 
SOLiD Total RNA-Seq kit according to the manufacturer's 
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instructions with the following modifications: (1) for RNA 
fragmentation, 300 ng of rRN A- depleted RNA were incu- 
bated at 95°C for 5 min. Bioanalyzer traces were used to 
confirm that the largest proportion of RNA fragments were 
around 200 bp. (2) Bioanalyzer traces following reverse 
transcription were used to confirm that hybridization and 
ligation of RNA adapters and subsequent reverse transcrip- 
tion of the RNA to cDNA was successful. (3) Amplification 
of cDNA was carried out following reverse transcription, 
but preceding the size-selection step. Samples were barcod- 
ed by replacing the SOLiD 3' primer with the appropriate 
barcoding primer as directed in the SOLiD instructions. 
Size selection of DNA amplicons ranging from 200 to 
300 bp was performed using an E-Gel iBase system (SYBR 
Safe 2% SizeSelect gel) (Life Technologies, Carlsbad, CA). 
This step was repeated to ensure that adapter-adapter 
dimers were fully eliminated before sequencing. Paired-end 
sequencing was carried out by the University of Texas 
Genomic Sequencing and Analysis Facility on a Life 
Technologies SOLiD 5500x1 sequencing system. Targeted 
sequencing depth was six-million paired-end reads per 
sample. 

Read mapping and visualization and 
transcript quantification 

Filtering and alignment of the sequencing data was per- 
formed at the UTGSAF using the AB SOLiD BioScope 
Whole Transcriptome pipeline (vl.3), for whole-transcrip- 
tome RNA-seq analysis. Based on gene annotation infor- 
mation, the pipeline produced genome alignment results 
as a compressed binary version of the Sequence Alignment 
Map (BAM files). Mapped reads were visualized using 
BamView in Artemis 13.2.0 (Rutherford et al. 2000). 

RNA-seq data analysis 

To determine the transcriptional abundance for each 
gene, the number of reads that mapped within each 
annotated sequence was determined. The number of reads 
per kb of transcript per million mapped reads (RPKM) 
was used to normalize the raw data (Mortazavi et al. 
2008), and mean RPKM values were determined for both 
the WT and gacA mutant samples. The complete dataset 
from this study has been deposited at the National Center 
for Biotechnology Information (NCBI) Sequence Read 
Archive (SRA) with the Accession No. GSE43083. A ratio 
of the mean RPKM values (gacA mutant/WT) was com- 
puted for each gene and displayed by gene order from the 
chromosome origin of replication (Fig. 1). Comparisons 
were performed using a modified f-test and genes with 
differences in gene expression based on the ratio were 
considered for further analysis when the P-value was less 



than 0.05 (except as otherwise mentioned) and the 
expression ratio was > 2.0 or < 0.5 (Table S3) (Wang 
et al. 2011). 

qPCR methods and analysis 

qPCR was performed at the Texas A&M GTL using a pre- 
viously described method (Wang et al. 2012a,b). RNA 
was reverse transcribed using random primers (Invitro- 
gen) and Superscript III (Invitrogen) at 50°C for 1 h and 
inactivated at 75°C for 15 min. SYBR Green reactions 
were performed using the ABI 7900 HT Fast System 
(Applied Biosystems, Foster City, CA) in 384 well optical 
reaction plates. Aliquots (1 /.(L) of cDNA (2 ng/reaction) 
or water (no-template control) were used as template for 
qPCR reactions with Fast SYBR Green PCR Master Mix 
(Applied Biosystems) and primers (500 nmol/L 
final concentration). Primer pairs gacSRTl-gacSRT2, 
gacART 1 - gacART2 , rpeBRTl-rpeBRT2, pipRTl-pipRT2, 
phzYRTl-phzYRT2, phzRRT 1 -phzRRT2 , phzIRTl- 
phzIRT2, hfqRTl-hfqRT2, rpoDRTl-rpoDRT2, iopBRTl- 
iopBRT2, 16SRT1-16SRT2, and rpoSRTl-rpoSRT2 were 
used to detect the expression of gacS, gacA, rpeB, pip, 
phzY, phzR, phzl, hfq, rpoD, iopB, 16S rRNA, and rpoS 
genes, respectively (Table S2). qPCR amplifications were 
carried out at 50°C for 2 min, 95°C for 10 min, followed 
by 40 cycles of 95°C for 15 sec and 60°C for 1 min, and 
a final dissociation curve analysis step from 65 to 95°C. 
Two technical replicates of each of three biological repli- 
cates were used for each experiment. Gene expression lev- 
els were analyzed using the relative quantification (AACt) 
method (Wang et al. 2011, 2012a,b). The housekeeping 
sigma factor rpoD was used as the reference gene to nor- 
malize our samples (ACt = Ct target — Ct rpoD ). A relative 
quantification (RQ) value was calculated as 2 exp- 
(AACt = ACt, arget — ACt reference ) for each gene. A P-value 
was computed using a f-test to measure the significance 
associated with each RQ value. Variations were consid- 
ered statistically significant when the P-value was <0.05. 
RQ values for the gacA mutant were then normalized to 
those of WT. 

Cloning of rsmA, rsmE, rsmX, rsmY, rsmZ, 
and phzR genes 

The rsmA and rsmE flanking sequences were used to 
design primers (rsmAl-rsmA2 and rsmEl-rsmE2) to 
amplify these genes and their promoter sequences. Fol- 
lowing amplification, DNA fragments and vectors 
(pLAFR3) were digested by £coRI and BamHl and ligated. 
The resulting plasmids were verified by sequencing and 
designated as pLAFR3-rsmA and pLAFR3-rsm£, respec- 
tively. Plasmids were introduced into P. chlororaphis 
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Pchl3084 3493-*-* 

tRNA-Ala -y* 




Gene order 



Figure 1. Differential gene expression between wild-type (WT) Pseudomonas chlororaphis 30-84 and a gacA mutant. Each point represents one 
of the 5980 annotated genes in the 30-84 genome, with the x-axis showing gene order (from the DNA replication origin), and the y-axis showing 
the log 2 of transcript abundance for each gene in the gacA mutant relative to the WT strain. The arrows point to a few of the genes or gene 
clusters that are differentially expressed. Pchl3084_3493, conserved hypothetical protein; tRNA-Ala, alanine transfer RNA; T6SS, type VI secretion 
system; rpoS, stationary phase sigma factor; rsmX, rsmY, rsmZ: ncRNAs; AprX, metallopeptidase; AprA, metalloprotease. 



strains by triparental mating as described previously 
(Wang et al. 2012b). Transformants were selected on LB 
plates supplemented with appropriate antibiotics. 

Primers pairs rsmXl-rsmX2, rsmYl-rsmY2, rsmZl- 
rsmZ2, and phzRl-phzR2 were used to amply rsmX, 
rsmY, rsmZ, and phzR genes. Primers rsmXl, rsmYl, 
rsmZl, and phzRl contain constitutive-hybrid P tac pro- 
moters (De Boer et al. 1983) that ensure the correct tran- 
scription start site. The amplified DNA fragments and 
vectors (pLAFR3 for ncRNAs and pPROBE-GT2 for 
phzR) were digested by EcoRl and BamHl and ligated 
resulting in plasmids pLAFR3-P MC -r5mX, pLAFR3- 
V tac -rsmY, pLAFR3-P tac -rsmZ, and pGT2-P tac -phzR. The 
genotypes were confirmed by both enzymatic digestion 
and sequencing. 

Construction of an rsmE mutant 

A 30-84 genomic library was introduced into strain 
30-84ZN (phzB::lacZ) via triparental mating. Exconjugants 
were screened for alterations in phenazine gene expression 
as indicated by differences in /?-galactosidase activity on 
medium containing X-gal. Cosmid pLSP298 caused a visi- 
ble reduction in phzB expression and was chosen for fur- 
ther analysis. Cosmid pLSP298 was mutated using EZ::TN 
<KAN-2> (Epicentre Biotechnologies, Madison, WI) and 
cosmids containing insertions were introduced into strain 
30-84ZN via triparental mating and compared with strain 
30-84ZN (pLSP298) for phzBr.lacZ expression. One EZ:: 
TN insertion (pLSP298EZ::TN5-ll-2) that no longer 
reduced phzB expression was cloned as a BamHl fragment 
containing the EZ::TN Km R gene into pIC20H (Marsch 
et al. 1984). The regions adjacent to the cloned region of 
the EZ::TN were sequenced using the EZ::TN supplied 
primer and shown to be rsmE. Cosmid pLSP298EZ:: 
TN5-11-2 was introduced into strain 30-84ZW, a AgacA 



/5-galactosidase reporter (phzB::lacZ, gacA~) via triparental 
mating. As a result, the rsm£::EZ::TN insertion was intro- 
duced into the genome of 30-84ZW via marker exchange 
to create a AgacA /J-galactosidase reporter with a muta- 
tion in rsmE, designated 30-84ZWE (phzB::lacZ, gacA~, 
rsmE~). 

Quantification of rsmZ promoter activity 

A 630 bp Sail fragment immediately upstream of rsmZ 
that contained the rsmZ promoter (V rsmZ ) was cloned into 
the Sail site within pIC20H. A 681 bp EcoRl-Hindlll 
fragment from the resulting plasmid was cloned into 
pPROBE-KT2 upstream of the promoter-less GFP gene to 
create pKT2-F rsmZ -gfp. A Hindlll fragment containing a 
promoter-less lacZ gene was cloned from plasmid 
pKOK6.1 and inserted into the unique HmdIII site 
between V rsmZ and gfp in pKT2-P rsmZ -gfp resulting in 
pKT2-P rsmZ -ZacZ-g/p. Plasmid pKT2-P rsmZ -lacZ-gfp was 
introduced into the 30-84Ice strain and gac mutants via 
triparental mating. GFP expression for strains containing 
pV rsmZ -lacZ-gfp was determined visually and quantified by 
/J-galactosidase activity (Miller 1972). 

Quantification of phenazine production 

Pseudomonas chlororaphis strains were grown with aera- 
tion at 28°C in LB for 24 h. Phenazines were extracted 
and quantified by UV-visible light spectroscopy as 
described previously (Wang et al. 2012b). Briefly, tripli- 
cate 10 mL cultures grown overnight at 28°C were centri- 
fuged (5000g), and the supernatants were acidified to ca. 
pH 2 with concentrated HC1. Phenazines were extracted 
with an equal volume of benzene for 6 h. Following 
separation and evaporation of the benzene phase under a 
stream of air, phenazines were resuspended in 0.5 mL of 
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0.1 N NaOH, and serial dilutions were quantified via 
absorbance at 367 nm. The absorbance for each sample 
was normalized to the total absorbance of the 10 mL 
culture. Strain 30-84ZN, used as a negative control, 
carries a mutation in the phzB gene and is deficient in 
phenazine production. 

AHL extraction and biological assays 

Total AHL extractions were prepared from cell-free super- 
natants of strain 30-84, AgacS, or AgacA with plasmid 
pLAFR3 containing either no insert (NI) or rsmZ. The 
cell-free supernatants were made from overnight cultures 
grown in LB medium (OD 600 ~1.8) as described previously 
(Pierson et al. 1998). Briefly, 5 mL cultures were grown 
overnight at 28°C with shaking in LB broth to an OD 600 
of 1.8. The cultures were centrifuged and supernatants 
were mixed with an equal volume of acidified ethyl ace- 
tate. The ethyl acetate phase was evaporated and the dried 
extracts containing AHLs were suspended in a volume of 
LB equal to the original culture, which was subsequently 
filter-sterilized (0.45 /im). AHL production was quantified 
by inoculating the extracted AHLs with the AHL-specific 
reporter P. chlororaphis strain 30-84I/Z (phzF phzB:: 
lacZ). Strain 30-84I/Z is deficient in AHL production, but 
responds to exogenously added AHL's by producing 
/5-galactosidase. /?-galactosidase activity was determined 
subsequently on cultures grown shaking at 28°C after 
24 h (Wood and Pierson 1996). The assays were repeated 
at least three times. 

Phenotypic analyses and fungal inhibition 
assays 

The production of siderophore, exoprotease, gelatinase, 
amylase, lipase, and HCN was qualitatively measured 
(Chancey et al. 1999, 2002). Strains were spotted onto 
appropriate media: KMB plates + FeCl 3 -6H 2 0 (final con- 
centration 1 mmol/L) for fluorescence, skim milk plates 
for exoprotease, gelatin plates for gelatinase, starch plates 
for amylase, and tributyrin plates for lipase, respectively. 
For HCN production, a 1-inch piece of cyantesmo 
cyanide indicator paper (Macherey-Nagel, Germany) was 
attached to the inside lid of sealed plates. Cell growth and 
phenotypes were measured after 24 and 48 h. For patho- 
gen inhibition, overnight cultures of 30-84 and derivatives 
were spotted onto triplicate LB + 0.5% potato dextrose 
agar plates. After 2 days of growth at 28°C, a 5 mm plug 
of G. graminis var. tritici was placed in the center of the 
plates. After 4 days, zones of inhibition, the distance 
between the edge of the bacterial colony and the fungal 
mycelium, were measured (Whistler and Pierson 2006). 
The assays were repeated twice. 



Bacterial swimming motility assays 

For P. chlororaphis WT, gac mutants, and 30-84Z strain, 
bacterial cell suspensions were grown overnight in LB 
broth. Five microliter of the bacterial suspension was pla- 
ted onto the center of motility agar plates (10 g tryptone, 
5 g NaCl, 2.5 g agar per liter distilled water) as previously 
described (Wang 2011). Diameters were determined fol- 
lowing incubation at 28°C for 48 h. The experiments 
were repeated at least three times. 

Results and Discussion 

Transcriptomic analyses of phenazine 
biosynthetic and regulatory genes 

The RNA-seq analysis identified 771 genes that were differ- 
entially expressed in the gacA mutant compared with the 
WT (Fig. 1; Table S3). Among these genes, the transcripts 
of 551 genes were underrepresented and 220 genes were 
overrepresented in the gacA mutant relative to WT. Genes 
encoding components of the Gac signal transduction path- 
way were among those differentially regulated (Fig. 1). 

Transcript abundance patterns in P. chlororaphis 
30-84 WT 

To better understand the role of the Gac system on phen- 
azine regulation, transcript abundance of the phenazine 
biosynthetic and regulatory genes in the WT were mea- 
sured. The mean RPKM determined by RNA-seq for gacA 
and gacS in the P. chlororaphis 30-84 WT were 266.9 and 
29.2 (Table 1), respectively, demonstrating that in the 
WT, gacA expression is about ninefold higher than gacS 
under these experimental conditions (e.g., grown in 
AB + 2% CAA to OD 600 = 1.2). RPKM values for rsmA 
and rsmE in WT were similar to each other (151.8 and 
292.7, respectively), and to gacA, whereas values for the 
three ncRNAs rsmX, rsmY, and rsmZ were significantly 
higher. In P. protegens Pf-5, GacA transcriptionally acti- 
vates rsmX, rsmY, and rsmZ by binding to conserved 
motifs in their promoter regions (Humair et al. 2010). 
Analysis of rsmX, rsmY, and rsmZ in P. chlororaphis 30-84 
revealed the conserved palindromic consensus sequence 
typical of GacA-controlled ncRNA genes and a conserved 
— 10 hexamer (TAATCT) promoter sequence identical to 
that found in P. protegens Pf-5 (Humair et al. 2010) 
(Fig. 2A). Analogous to P. protegens Pf-5, the rsmX and 
rsmY promoter sequences are similar in length, but the 
rsmZ promoter sequence is longer (Fig. 2A). Consistent 
with rsmY having a promoter with the highest similarity 
to the conserved GacA binding site, rsmY was expressed 
at the highest level (RPKM = 38112.4), approximately 
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Table 1. Mean transcript abundance and ratio of abundances (Agac/VWT) of the phenazine biosynthetic and regulatory genes in the gacA 
mutant compared with the WT strain. 



Gene ID 




Gene 


Protein description 


Mean RPKM WT 


Mean RPKM AgacA 


Agac/VWT 


P-value 


gacS/gacA 
















Pchl3084_3491 


gacA 


Response regulator 


266.9 


1.43 


0.01 


0.00 


Pchl3084_4333 


gacS 


Sensor protein 


29.17 


21.15 


0.73 


0.05 


phz genes 
















Pchl3084_4951 


phzX 


Phenazine biosynthesis protein 


40.50 


4.52 


0.11 


0.03 


Pchl3084_4952 


phzY 


Phenazine biosynthesis protein 


37.87 


5.28 


0.14 


0.03 


Pchl3084_4953 


phzF 


Phenazine biosynthesis protein 


20.00 


2.63 


0.13 


0.03 


Pchl3084_4954 


phzA 


Phenazine biosynthesis protein 


8.52 


0.78 


0.09 


0.15 


Pchl3084_4955 


phzB 


Phenazine biosynthesis protein 


13.47 


2.08 


0.15 


0.04 


Pchl3084_49S6 


phzC 


Phenazine biosynthesis protein 


8.24 


1.90 


0.23 


0.00 


Pchl3084_49S7 


phzD 


Phenazine biosynthesis protein 


7.63 


1.40 


0.18 


0.10 


Pchl3084_49B8 


phzO 


Phenazine biosynthesis protein 


44.23 


21.70 


0.49 


0.06 


Regulatory genes 














Pchl3084. 


4949 


phzl 


Autoinducer synthase 


5.24 


1.50 


0.29 


0.01 


Pchl3084. 


4950 


phzR 


Transcriptional activator 


125.1 


67.12 


0.54 


0.02 


Pchl3084. 


.1189 


rpoS 


Sigma factor 


1045.7 


101.8 


0.10 


0.03 


Pchl3084. 


.51 BB 


pip 


Phenazine-inducing protein 


123.4 


105.4 


0.85 


0.04 


Pchl3084. 


.1659 


iopA 


Inducer of phenazine A 


49.10 


8.35 


0.17 


0.01 


Pchl3084. 


.1660 


iopB 


Inducer of phenazine B 


20.71 


2.51 


0.12 


0.00 


Pchl3084. 


.3224 


rpeA 


Sensor histidine kinase 


46.78 


49.28 


1.05 


0.67 


Pchl3084. 


.3225 


rp@B 


Cocrvinto ream i ^a+fur 
r\fcrS|JU[ lie [cyuldLUr 


45 37 


55 20 


1 22 


0 07 


Pchl3084. 


.2449 


csal 


Autoinducer synthase 


8.15 


8.72 


1.07 


0.65 


Pchl3084. 


.2451 


csaR 


Transcriptional activator 


0.54 


0.38 


0.71 


0.50 


Pchl3084. 


.0554 


hfq 


RNA chaperone 


1642.19 


1408.2 


0.84 


0.28 


Pchl3084. 


.4387 


rsmA 


Translational regulator 


151.82 


307.24 


2.02 


0.03 


Pchl3084. 


.2024 


rsmE 


Translational regulator 


292.7 


204.83 


0.70 


0.06 


Pchl3084. 


.3970 


rsmX 


ncRNA 


1486.51 


201.13 


0.14 


0.02 


Pchl3084. 


.5419 


rsmY 


ncRNA 


38112.4 


329.2 


0.01 


0.02 


Pchl3084. 


.1190 


rsmZ 


ncRNA 


1679.64 


261.26 


0.16 


0.05 



22-fold higher than rsmX (RPKM = 1486.5) or rsmZ 
(RPKM = 1679.6) in the WT (Table 1). 

We also looked at the relative transcript abundance of 
other regulatory genes associated with phenazine produc- 
tion or quorum sensing (Table 1). In the WT, hfq and 
rpoS were expressed at very high levels with RPKM values 
of 1642.2 and 1045.7 (the average RPKM value is 119.5), 
respectively. These two regulators were of interest because 
the RNA chaperone, Hfq, has been shown previously to 
control Jux^-type regulators (Wilf et al. 2012) and the 
sigma factor, RpoS, to control phzR/phzI through pip (Gir- 
ard et al. 2006a; Wang et al. 2012b). RPKM values for the 
two essential phenazine regulators pip and phzR were 
123.4 and 125.1, around 10-fold lower than those of hfq 
and rpoS. The RPKM value for the other quorum sensing 
regulator csaR, which plays a minor role in phenazine reg- 
ulation under laboratory conditions, was only 0.54 sug- 
gesting it was barely expressed under the assay conditions. 
The RPKM values for both quorum sensing signal synth- 
ases, phzl and csal, were 5.2 and 8.2, respectively. 

Similar to other phenazine-producing bacteria, the 
core phenazine, PCA, is encoded by a conserved set of 



biosynthetic genes (Mavrodi et al. 2006; Mentel et al. 
2009); in 30-84 these are phzXYFABCD (Pierson et al. 
1995). According to the RNA-seq analysis, the transcrip- 
tion start site of the phzX was predicted 115 nt upstream 
of the translation start site (ATG) (Fig. 2B). This was 
consistent with the expected structure of a quorum sens- 
ing-regulated promoter containing a —10 hexamer and a 
Lux-box located at the —35 site predicted for PhzR 
binding (Qin et al. 2007) (Fig. 2C). The transcript 
abundance of the phz genes showed a staircase-like 
pattern from phzX to phzD (Table 1), suggesting 
decreased transcription of genes with position along the 
phz operon. In many phenazine-producing Pseudomonas 
species, PCA is then modified into other derivatives by 
one or more accessory genes encoding modifying enzymes 
(Mavrodi et al. 2006). In 30-84, PCA is modified into the 
2-hydroxy forms, 20HPCA and 20HPZ, by an enzyme 
encoded by phzO, adjacent to the phenazine biosynthetic 
operon (Maddula et al. 2008). It was shown previously 
that the expression of phzO was driven by its own pro- 
moter (Maddula et al. 2008) and this was supported by 
the RNA-seq analysis (data not shown). The expression of 
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(A) 

Upstream Activation Sequence (UAS) 

ITGTAAG NN ATN NCTTACAl 



rsm X UAS 

GATC TATGGCGAAAGGCTTACAT GCGCTGGCTGCTATCTCGCCTGTTGCC 
GATTGGGCGACAGGC fTAATCt TCACC 

rsmY UAS 

gtc a 1 tgtaag c attctcttac a^ tg acgtg ctg ctaatcaccc atatcg g c a 

gattggatccagag BaatctA cttcac 
rsmZ 

UAS 

AGC A'TGTAAGTAAAGGCTGACTt TTTCGTAAGTGTTCGGTTTTTATGCGGT 
ATGAAAGTTGTCTGTTTGATCTGAGGTTTTTAGTTATCTTATTGATTAATAA 
GAACTTATTTTAATGATTGGGGCGGTGTTCGACACGTTTTCTACCTGTCAG 
TGCCATTGCTCTGTCATGGAAACTCTC TAATATCt AGCCCTG 




CCACTC 



ATG 



^ 35 pfebOX 10 

aatc Kctacaagatctggtagi1 tccacccccaagaaatgcaggtg |tataaa| 
+1 

caacacccactcagcaccgccacgcaggaacaatagttaaacgctgtaaa 



tattcacccatgtcagttaactaaatatttactttcaacactatcaccccca 
acta aggagg atgctgcccatgcctgcttcgctttcccctagc 

Figure 2. Promoter regions of ncRNA and phenazine biosynthetic 
genes. (A) Alignment of the rsmX, rsmY, and rsmZ promoter regions of 
Pseudomonas chlororaphis 30-84. The consensus Upstream Activation 
Sequence recognized by GacA (UAS) (Humair et al., 2010) is shown. The 
putative UAS for each promoter region is boxed and conserved 
nucleotides are in bold. The -10 promoter elements are also indicated. 
(B) RNA-seq profile showing sequence reads across the phz promoter 
region. RNA-seq results were visualized using the Artemis 
(http://www.sanger.ac.uk/resources/software/artemis/) genome browser. 
Transcription and translation start sites are shown below the profile. (C) 
Nucleotide sequence of the phz promoter. Boxed nucleotides indicate 
sequences involved in transcription initiation (PhzR Box, RNA polymerase 
binding site), whereas sequences involved in translation initiation 
(ribosome or RsmE binding sites) regulatory elements are underlined. 
Arrow indicates the direction of transcription from the start site. 



phzO was similar to that of phzX (RPKM values of 40.5 
and 44.2, respectively, Table 1). 

Phenazine biosynthetic and regulatory gene 
transcripts are in lower abundance in the gacA 
mutant relative to WT 

In order to identify genes regulated by GacA in 30-84, 
mean RPKM values for the gacA mutant were compared 
with WT values. As expected, the RPKM value for gacA 



was very low in the gacA mutant compared with the WT 
(1.43 vs. 266.9, respectively) demonstrating that the gac A 
message was successfully disrupted, whereas expression of 
gacS was unchanged (Table 1). Mutation in gacA also 
resulted in lower abundance of the three ncRNAs; the 
greatest fold change was for rsmY (~1 15-fold decrease). 
The RPKM values for rsmX, rsmY, and rsmZ were 
decreased to 201.1, 329.2, and 261.3, respectively 
(Table 1). These data indicate that GacA positively regu- 
lates the expression of rsmX, rsmY, and rsmZ in P. chloro- 
raphis 30-84. Interestingly, the transcript abundance of 
rsmA was 2.2-fold higher in the gacA mutant compared 
with WT 30-84, whereas that of the closely related rsmE 
was relatively unchanged (Table 1). 

Consistent with the deficiency in phenazine produc- 
tion observed in gac mutants, mutation in gacA resulted 
in 2- to 11 -fold decreases in the transcript abundance of 
the phenazine biosynthetic genes (Table 1). Transcripts 
of six regulatory genes involved in phenazine produc- 
tion, iopA, iopB, rpoS, pip, phzR, and phzl, were under- 
represented by one- to ninefold in the gacA mutant as 
compared with WT, whereas other regulatory genes such 
as hfq, csal/R, and rpeA/rpeB were not appreciably 
altered (Table 1). The expression of these phenazine reg- 
ulatory and biosynthetic genes was verified by qPCR 
(Fig. 3). The results suggest that GacA controls the 
expression of the phenazine biosynthesis genes in 30-84 
partly by activating the expression of phzR, phzl, pip, 
rpoS, and iopAliopB. 




gac A gacS rpoS iopB rpeB pip phzR phzl phzY 



Figure 3. Verification of RNA-seq data by quantitative reverse 
transcription-polymerase chase reaction (qPCR). The relative fold 
change of gacA, gacS, rpoS, iopB, rpeB, pip, phzl, phzR, and phzB 
genes in gacA mutant compared towith the WT strain. The gene 
rpoD was used as the reference. Cells were grown in AB minimal 
medium + 2% casamino acid for 18 h with shaking to an OD 60 o of 
1.2. Data points represent means of three replicates ± standard 
deviations. These experiments were repeated at least three times and 
similar results were obtained. 
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Genetic analyses of the phenazine 
biosynthetic and regulatory genes in 30-84 

Expression of phzR/phzl does not rescue 
phenazine production in a gacA mutant 

It was shown previously that the expression of phzl and 
AHL production are significantly reduced in 30-84 gacS 
and gacA mutants (Chancey et al. 1999). RNA-seq analy- 
sis indicated that both quorum sensing genes phzl and 
phzR were significantly decreased in the gacA mutant. 
One hypothesis is that this deficiency in quorum sensing 
is the reason for the lack of phenazine production by 
gacA mutants. However, if Gac operates by controlling 
phenazine production posttranscriptionally, constitutive 
expression of quorum sensing genes should not comple- 
ment a gac mutant. To test this, the phzR gene was cloned 
into a medium copy vector driven by the V tac promoter, 
and the plasmid with or without the phzR insertion was 
introduced into either WT, 30-84R (phzR~), or the gacA 
mutant. As expected, complementation of the phzR 
mutant by the pGT2-P tac -phzR plasmid or the phzl 
mutant 30-841 by the addition of AHLs fully rescued 
phenazine production (Fig. 4A and B). In contrast, 



(A) 




WT AphzR AphzR Aphzl Aphzl AgacA 

(NI) (NI) (phzR) (-AHL) (+AHL) (phzR+ARL) 



(B)4 




WT AphzR AphzR Aphzl Aphzl AgacA 

(NI) (NI) (pfeff) (-AHL)(+AHL)(/jte^+AHL) 



Figure 4. Phenazine production by WT and mutants with additional 
copies of phzR and supplemental AHLs. (A) Bacterial strains were 
grown in 2 ml_ LB medium in a 24-well plate at 28°C without shaking. 
The photo was taken 48 h after inoculation. (B) Bacterial strains were 
grown in LB medium for 24 h at 28°C with shaking at 200 rpm. 
Treatments included strains with the plasmid pGT-2 containing either 
no insert (NI) or the plasmid containing a ? tac -phzR insertion, AHLs 
were added to the growth medium as indicated. Phenazine was 
extracted and quantified at OD 367 . AHLs were extracted as in ethyl 
acetate as described previously and total AHL quantified by ability to 
rescue phenazine production in the phzl mutant 30-841. Data points 
represent means of three replicates ± standard deviations. 



constitutive expression of phzR with the addition of AHLs 
to the growth medium did not rescue phenazine produc- 
tion in the gacA mutant (Fig. 4 A and B). These results 
demonstrate that the underexpression of the PhzR/PhzI 
quorum sensing system in Gac mutants is not solely 
responsible for the phenazine deficiency observed. 

RsmE, but not RsmA, regulates phenazine and 
AHL production 

Pseudomonas chlororaphis 30-84 has both of the RNA-binding 
proteins RsmA and RsmE (this study) and they share 72% 
and 87% sequence identity and similarity, respectively. Amino 
acid identity of 30-84 RsmA and RsmE to their respective 
homologs in other Pseudomonas species was high (e.g., 100% 
and 93% identical to P. protegens Pf-5, P. fluorescens PfOl, 
P. syringae DC3000, P. aeruginosa PAOl; 94% and 91% iden- 
tical Escherichia coli 0157). The current model predicts that 
overexpression of RsmA or RsmE in the WT should provide 
phenotypes that are similar to gac mutants (Heeb and Haas 
2001). To determine the roles of RsmA and RsmE in the reg- 
ulation of phenazine production, multiple copies of rsmA and 
rsmE were introduced separately into the WT strain. Phen- 
azine production decreased 7.1 -fold in the WT strain harbor- 
ing pLAFR3-rsm£ (compared to WT with pLAFR3-NI) 
whereas overexpression of rsmA did not alter phenazine 
production in the WT strain (Fig. 5A). To confirm that rsmA 
was overexpressed in the WT, qPCR was performed and 
results showed that multiple copy of rsmA leads to 
12.8 ± 0.6-fold transcript increase in the WT strain. To cor- 
relate phenazine production with phenazine gene expression, 
phzBdacZ expression was measured in 30-84ZN harboring 
these plasmids. Consistent with phenazine production, multi- 
ple copies of rsmE reduced 30-84ZN /J-galactosidase activity 
from 487.0 ± 78.6 to 48.1 ± 5.2 in 30-84ZN with 
pLAFR3-NI compared with pLAFR3-rsm£, respectively; the - 
galactosidase activity produced by strain 30-84ZN harboring 
pLAFR3-rsmA was similar to the 30-84ZN (pLAFR3-NI) con- 
trol (511.0 ± 42.1 vs. 487.0 ± 78.6, respectively) (Fig. 5B). 
The result that RsmE, but not RsmA is involved in the regula- 
tion of phenazines is consistent with the observation in Pseu- 
domonas sp. M18 that RsmA negatively controls pyoluteorin 
production but not phenazine production (Zhang et al. 
2005). As pyoluteorin genes were not detected in P. chlorora- 
phis 30-84, the target(s) of RsmA remains to be discovered. 

To determine whether RsmE also inhibited AHL 
production, the effect of overexpression of rsmE on AHL 
signal production was measured. Total AHL signal was 
extracted from 24 h bacterial cultures of WT (pLAFR3- 
NI) and WT (pLAFR3-rsm£). The relative amount of AHL 
signal was quantified by activation of /^-galactosidase 
activity in strain 30-84IZ (Fig. 5C). /3-galactosidase activity 
of strain 30-84IZ was approximately eightfold lower when 
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Figure 5. Role of RsmE in AHL and phenazine production. (A) Phenazine production by 30-84 with plasmid pLAFR3 containing either no insert 
(NI), rsmA, or rsmE. Bacterial strains were grown in LB medium for 24 h at 28°C with shaking. Phenazine was extracted as described previously 
and total phenazine was quantified at OD 367 . Data points represent means of three replicates ± standard deviations. (B) Repression of phzB 
expression in 30-84ZN by introduction of extra copies of rsmE but not by introduction of extra copies of rsmA. Treatments included 30-84ZN with 
plasmid pLAFR3 containing either NI, rsmE, or rsmA. 30-84ZN is a translational phzB::lacZ fusion. The /?-galactosidase activities (Miller Units) were 
determined in triplicate (mean ± SE). (C) Suppression of AHL production by rsmE in trans. AHL accumulation was quantified from overnight 
cultures of WT with pLAFR3 containing either NI or rsmE using the AHL-specific reporter 30-841/Z (phzl~, phzB::lacZ). The relative amount of AHL 
was determined by /?-galactosidase assays. (D) Effect of inactivation of rsmE on phzB expression. Treatments included the phzB::lacZ reporter (30- 
84ZN), a phzB::lacZ/GacA" derivative (30-84ZW), and a p/7zfi.\7acZ/GacA~/RsrnE _ derivative (30-84ZWE). Bacterial cultures were grown in LB and 
phzB expression was measured as /?-galactosidase activity after 24 h. Data are the means of two replicates from one representative experiment. 



the reporter strain was grown in the presence of AHLs 
prepared from overnight cultures of WT (pLAFR3-rsm£) 
compared with AHLs prepared from WT (pLAFR3-NI). 
These results indicate that RsmE overexpression also 
resulted in a significant decrease in AHL production. 

To determine whether RsmE functions downstream of 
GacA in the regulation of phenazine gene expression, 
cosmid pLSP298EZ:TN5-ll-2 containing a transposon 
insertion that inactivated rsmE was introduced into the 
genome of the gacA mutant 30-84ZW via marker 
exchange mutagenesis. Comparison of phzB::lacZ 
expression in strains 30-84ZN, 30-84ZW, and 30-84ZWE 
indicated that loss of rsmE restored /3-galactosidase 
activity to the gacA mutant (Fig. 5D). 

Constitutive expression of rsmZ rescues 
phenazine and AHL production in gac mutants 

We further examined whether constitutive expression of 
ncRNAs would remove the need for a functional Gac 



system for phenazine production in strain 30-84. The 
rsmX, rsmY, and rsmZ genes were cloned under the control 
of a constitutive P tac promoter such that the +1 site of the 
promoter was at the beginning of these genes. Although 
strains AgacS (pLAFR3-NI) and AgacA (pLAFR3-NI) failed 
to synthesize measurable amounts of phenazines, strain 
AgacS (pLAFR3 -P tac - rsmZ) and AgacA (pLAFR3-P tac -«mZ) 
produced approximately the same amount of phenazines 
as WT (pLAFR3-P fflc - rsmZ) (Fig. 6A). In contrast, constit- 
utive expression of rsmX did not rescue phenazine produc- 
tion (data not shown). The expression rsmX was 
confirmed using qPCR, the gacA mutant containing rsmX 
plasmid showed 2.4-fold increase of rsmX expression than 
the WT strain. Constitutive expression of rsmY by the P I£JC 
promoter appeared to be lethal in gac mutants, as no 
transformants were recovered with electroporation or 
triparental mating (data not shown). 

We further explored the regulatory role of GacA on 
ncRNA expression by characterizing the effect of 
mutations in gacA on rsmZ promoter activity using the 
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Figure 6. Impact of rsmZ on phenazine and AHL production. (A) Phenazine production by different derivatives in vitro. Treatments included 
strains with plasmid pLAFR3 containing either no insert (NI) or P ac -rsmZ. Growth conditions were as in (Fig. 5A). Phenazine was extracted in 
benzene and the amount of phenazine was measured at OD 367 . Data points represent means of three replicates ± standard deviations. Similar 
results were obtained in at least three independent experiments. (B) Promoter activity of the rsmZ gene in WT, gacS, and gacA was determined 
by /?-galactosidase activities. Treatments included strains with plasmid pKT-2 containing either a promoter-less lacZ insertion or the rsmZ promoter 
fused to the promoter-less lacZ gene (? rsmr lacZ). Insertions were made in front of the promoter-less gfp gene, already present on the plasmid. 
Bacterial strains were grown in LB medium for 24 h at 28°C with shaking to an OD 60 o of -1.8. Each value is the average from two different 
cultures ± standard deviations with two independent experiments. The experiment was repeated at least twice and similar results were obtained. 
(C) Effect of rsmZ on AHL production. Total AHL extractions were prepared from cell-free supernatants of WT, AgacS, or AgacA with plasmid 
pLAFR3 containing either NI or P tac -rsmZ. Cell-free supernatants were made from overnight cultures grown in LB medium (OD 60 o -1.8). The 
relative amount of AHL in each extract was determined by /?-galactosidase assays using the AHL-specific reporter 30-84I/Z (phzl~, phzB::lacZ). Bars 
represent /J-galactosidase activity measured in Miller Units of 30-84I/Z grown with the AHL extracts indicated. Each value is the average from 
three different cultures ± standard deviations. (D) Relative expression of phzl, phzR, and phzY genes in the gacA mutant (harboring the pLAFR3- 
? tac -rsrnZ) compared with WT (harboring pLAFR3 empty vector) as determined by qPCR. rpoD was used as the reference gene. Cells were grown 
in AB minimal medium + 2% casamino acid for 18 h with shaking to an OD 60 o of 1.2. Data points represent means of three 
replicates ± standard deviations. These experiments were repeated at least three times and similar results were obtained. 



reporter plasmid pKT2-P rsmZ -lacZ-gfp (e.g., by measuring 
fluorescence or /J-galactosidase activity). The reporter 
plasmid was introduced into a gacA and a gacS mutant as 
well as 30-84Ice (phzB::inaZ), which has a functional 
gacA, but does not produce phenazines (thus no inter- 
ference with GFP fluorescence). Strain 30-84Ice (pKT2- 
PrsmZ-lacZ-gfp) was brightly fluorescent whereas the Gac 
mutants with the reporter plasmid were not (data not 
shown). Introduction of the plasmid into the gacA or 
gacS mutant resulted in an 11 -fold reduction in expres- 
sion of /J-galactosidase relative to strain 30-84Ice (pKT2- 
PrsmZ-lacZ-gfp) (Fig. 6B). These results confirm that the 
Gac system positively regulates the expression of rsmZ. 

To determine whether constitutive expression of rsmZ 
also regulated AHL signal production, the effect of 



constitutive expression of rsmZ on AHL signal production 
was measured. Total AHL signal was extracted from 24 h 
bacterial cultures and the relative amount of AHL signal 
present was quantified by activation of |8-galactosidase 
activity in strain 30-84IZ (Fig. 6C). The /3-galactosidase 
activity of the 30-84IZ reporter was significantly higher 
when grown with AHLs prepared from strains /AgacS 
(pLAFR3-P MC -rsmZ) and AgacA (pLAFR3-P fflc -rsmZ) com- 
pared with AHLs prepared from AgacS (pLAFR3-NI) and 
AgacA (pLAFR3-NI); /?-galactosidase activities of the 
30-84IZ reporter were 712.0 ± 77.1 and 723.1 ± 56.2 as 
compared with 28.9 ± 12.2 and 24.1 ± 10.2, respectively. 
This increase (27-fold) in phzB expression with constitu- 
tive rsmZ expression supports the hypothesis that the 
presence of the rsmZ ncRNA bypasses the requirement 
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for a functional Gac system for both phenazine and AHL 
production. qPCR analyses confirmed that the transcript 
abundance of phzl, phzR, and phzY in the gacA mutant 
was two- to fourfold higher in the presence of rsmZ 
(Fig. 6D). Collectively, these results indicate that Gac reg- 
ulates phenazine production and quorum sensing by acti- 
vating the expression of rsmZ. 

Constitutive expression of rsmZ also qualitatively restored 
HCN, exoprotease, gelatinase, and lipase production and 
reduced fluorescence to WT levels in strain 30-84W (data 
not shown). Consistent with restoration of these capabilities, 
constitutive expression of rsmZ also reinstated 30-84's ability 
to inhibit the mycelial growth of the fungal pathogen Ggt. 
Zones of Ggt mycelial inhibition by WT, AgacA, AgacA 
(pLAFR3-NI), and AgacA (pLAFR3-rsmZ) were 8.4 ± 0.9, 
0.0 ± 0.3, 0.2 ± 0.4, and 9.4 ± 0.8 mm, respectively. 

Model of Gac/Rsm regulation of phenazines 
in 30-84 

In this study, the major components of the Gac/Rsm signal 
transduction pathway in P. chlororaphis 30-84 were identi- 
fied and characterized (Fig. 7). For example, we discovered 
that the genome of 30-84 contains genes encoding both 
RsmE and RsmA. Only overexpression of rsmE in the WT 
reduced the gene expression and production of phenazine. 
Similarly only deletion of rsmE in a gacA mutant restored 
both, demonstrating that rsmE deletion bypassed the need 
for a functional GacA to produce phenazines. The genome 
of 30-84 contains the gene sequences for the three small 
ncRNAs, rsmX, rsmY, and rsmZ. A functional Gac system is 
required for rsmZ expression and that expression of rsmZ 
from a constitutive promoter restored phenazine, HCN, 
exoprotease, lipase, and gelatinase production and reduced 
fluorescence to WT levels in a gacA mutant. In contrast, 
constitutive expression of rsmX did not restore any of the 
pheno types tested, and constitutive expression of nor rsmY 
could not be achieved. The importance of the rsmZ in Gac 
regulation of products essential for biological control is 
interesting given that rsmZ has the longest promoter region, 
suggesting the potential for multiple protein interactions in 
its regulation. The observation that rsmZ plays a different 
role than the other two ncRNAs in Gac-mediated regulation 
has not been reported previously. However, the observation 
that transcript abundances of both rsmX and particularly 
rsmY are significantly lower in a GacA mutant is consistent 
with the hypothesis that they are involved in Gac-mediated 
regulation of as yet unidentified targets. Consistent with 
Gac/Rsm exerting posttranscriptional control of phenazine 
production, previous work showed that the introduction of 
a cosmid containing the entire phenazine biosynthetic clus- 
ter and phzR and phzl in trans on a low copy plasmid failed 
to complement a gacA mutant for phenazine production 
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Figure 7. Proposed model for the regulation of phenazine 
biosynthesis by the Gac-Rsm system in Pseudomonas chlororaphis 
strain 30-84. The model shows the presumed position of the GacS/ 
GacA system in relation to other known regulators of phenazine 
production. GacS is a transmembrane protein according to in silico 
predictions based on amino acid sequences. Amino acids H (histidine) 
and D (aspartate) involved in phosphorylation are indicated. The 
shaded area of GacA indicates the DNA-binding domain. Solid arrows 
and blunt lines point to genes (or processes) that are positively or 
negatively affected, respectively. Dashed arrows or lines indicate 
unknown or as yet uncharacterized regulatory pathways. GacA 
positively controls the expression of rsmX, rsmY, and rsmZ, which in 
turn activates phenazine production by titrating the translation 
suppressor RsmE. Although the transcript abundance of rsmA was 
2.2-fold higher in the gacA mutant compared with WT 30-84, RsmA 
does not appear to play a role in the regulation of phenazines under 
the experimental conditions tested here. A functional GacA is also 
required for increased expression of other phenazine regulatory 
genes, including rpoS, pip, and phzRIphzl. HK, histidine kinase; RR, 
response regulator; P, phosphoryl group. 



(Chancey et al. 1999). In this study, we bypassed transcrip- 
tional regulation of quorum sensing by providing PhzR 
(from a constitutive promoter) and purified AHL signal, 
confirming the prediction that Gac/Rsm controls phenazine 
gene expression posttranscriptionally. 

The transcriptomic data demonstrated that GacA controls 
the expression of the phenazine biosynthetic genes (phzXY- 
FABCD and phzO), the quorum sensing genes phzR and phzl, 
and also several additional regulatory genes involved in phen- 
azine gene expression including pip, rpoS, and iopA/iopB. In 
P. chlororaphis stains 30-84 and PCL1391, Pip promotes 
phenazine production by enhancing phzl and phzR expres- 
sion (Girard et al. 2006a; Wang et al. 2012b). In both strains, 
the expression of pip is regulated by the sigma factor rpoS, 
which in turn is regulated by GacA. In 30-84 pip also is regu- 
lated by the RpeB/RpeA TCST system (Wang et al. 2012b). 
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Complementation assays showed that constitutive expression 
of pip or phzR rescued phenazine production in rpeB 
mutants whereas multiple copies of rpeB or pip failed to res- 
cue phenazine production in a phzR mutant; for both strains, 
phzR expression restored phenazine production in a pip 
mutant (Wang et al. 2012b). Furthermore, expression of 
gacS, gacA, and rpoS were unchanged in an rpeB mutant, and 
conversely the expression of rpeA and rpeB were unchanged 
in an rpoS mutant (Wang et al. 2012b). These results indicate 
that RpoS and the RpeB/RpeA independently regulate phen- 
azine production using pip as a common regulatory interme- 
diate. Consistent with these findings, data presented here 
showed that expression of rpeB and rpeA were not altered by 
mutation of gacA, indicating that RpeB/RpeA and GacS/ 
GacA also do not regulate phenazines in a hierarchic manner. 
However, the effect of gacA mutation on rpoS and pip expres- 
sion suggests that GacA influences the expression of addi- 
tional phenotypes via RpoS and Pip. 

Transcriptomic analyses of genes controlled 
by Gac: a preliminary survey 

GacS/GacA functions as a positive regulator of 
secondary metabolism and exoenzyme production 

In 30-84, gac mutants are not only deficient in phenazine 
production, but also HCN, exoprotease, and lipase (Chancey 
et al. 1999). Consistent with these observations, the tran- 
script abundances of genes that encode for HCN, exopro- 
tease/exopeptidase, and lipase activities were lower in the 
gacA mutant than in WT (Table 2). Additionally, the tran- 
script abundances of 30-84 genes annotated to encode pyr- 
rolnitrin biosynthesis (prnABCD) and four genes involved in 
chitin degradation were lower in the gacA mutant than in 
WT (Table 2). Some of these results were confirmed by 
qPCR analysis using the hcnA, aprA, and prnA genes (data 
not shown). These results are consistent with previous find- 
ings showing a reduction in the abundance of transcripts of 
the biosynthetic genes for HCN, pyrrolnitrin, and protease 
in a P. protegens Pf-5 gacA mutant (Hassan et al. 2010). In 
Pf-5, the products of these genes have been demonstrated to 
be important for biological control activity. 

Transcripts of genes potentially involved in type 
VI secretion system, signal transduction, and 
carbon storage are in lower abundance in the 
gacA mutant relative to WT 

Genes annotated as encoding components of a type VI 
secretion system (T6SS) also were notably underexpressed 
(2- to ~100-fold) in the gacA mutant relative to WT, sug- 
gesting positive control by GacA (Table S4). Again, this is 
similar to previous findings showing that the expression of 
Pf-5 genes homologous to the P. aeruginosa T6SS (HIS-I 



type) was reduced in a gacA mutant from early to late 
growth phase (Hassan et al. 2010). In contrast to Pf-5 
where all 16 of the differentially expressed genes encoding 
the T6SS system were located in the same gene cluster, the 
24 differentially expressed T6SS genes identified in this 
study were distributed among six different loci. Previous 
studies also demonstrated GacS/GacA control of the T6SS 
in P. aeruginosa and P. syringae (Silverman et al. 2012). In 
P. aeruginosa, the T6SS delivers bacteriolytic effectors to 
target cells and is associated with competition among bac- 
terial species (Russell et al. 2011; Julie et al. 2012). Future 
studies will be required to determine whether 30-84 
contains a functional T6SS and if so, its function. 

Genes annotated as encoding enzymes involved in poly- 
hydroxyalkanoate (PHA) synthesis and utilization also 
were underexpressed (2- to ~100-fold) in the gacA mutant 
relative to WT, suggesting positive control by GacA 
(Table S5). PHAs are bacterial storage materials produced 
by many fluorescent pseudomonads including P. aerugin- 
osa, P. putida, and P. chlororaphis (aureofaciens) under 
nutrient conditions in which carbon is in relatively high 
abundance, but other nutrients are limited (Nishikawa 
et al. 2002; Pham et al. 2004; Hervas et al. 2008). Genes 
included phaD, a TetR-like transcriptional regulator that 
in P. putida KT2442 activates the pha cluster 
(phaClZC2D) conferring the ability to degrade PHAs (de 
Eugenio et al. 2010). In P. oleovorans GPol, phaD is tran- 
scribed as part of the pha operon and controls its own 
transcription as well as the transcription of the phalF 
operon, which in turn is involved in the regulation of pha 
genes (Prieto et al. 1999). Consistently, in this study, the 
transcripts of phaCl, phaZ, phaC2, phal, and phaF also 
were underrepresented (2- to 17-fold) in the gacA mutant 
relative to the WT (Table S5). Previously, it was shown 
that the PHA synthase phaC2 is underexpressed in an 
rpoS mutant of P. chlororaphis PCL1391 (Girard et al. 
2006b). Work with P. aeruginosa PAOl demonstrated that 
mutants deficient in PHA biosynthesis were less tolerant 
of heat shock both in liquid culture and in biofilms and 
were altered in attachment and biofilm architecture rela- 
tive to WT (Pham et al. 2004). Moreover, comparison of 
alginate production by PHA-negative derivatives of WT 
PAOl and the alginate-overproducing strain, FRD1, led 
the authors to suggest that the alginate and PHA path- 
ways compete for acyl-CoA precursors (Pham et al. 
2004). In 30-84, the transcript abundance of the alginate 
regulatory gene algK was fourfold higher in the gacA 
mutant compared with WT, however, the majority of 
genes annotated as being involved in alginate synthesis 
were unchanged. The roles of PHA or alginate production 
and utilization by 30-84 have yet to be determined, but it 
is interesting to consider given PHA's potential for influ- 
encing bacterial survival and biofilm formation. 
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Table 2. Mean transcript abundance and ratio of abundances (AgacAANT) of genes involved in secondary metabolism and iron uptake in the 
gacA mutant compared with WT strain. 

Gene ID Gene Protein description Mean RPKM WT Mean RPKM AgacA AgacAANJ P-value 



Hydrogen cyanide 



Pchl3084. 


.2378 


hcnA 


Hydrogen cyanide synthase 


746.54 


243.92 


0.33 


0.01 


Pchl3084. 


.2379 


hcnB 


Hydrogen cyanide synthase 


474.43 


70.53 


0.15 


0.00 


Pchl3084. 


.2380 


hcnC 


Hydrogen cyanide synthase 


705.07 


66.39 


0.09 


0.00 


Protease 
















Pchl3084. 


.3127 


aprA 


Metalloprotease 


100.83 


4.77 


0.04 


0.02 


Pchl3084. 


.2293 


aprX 


Metallopeptidase 


164.03 


6.44 


0.04 


0.01 


Lipase 
















Pchl3084. 


.3120 




Lipase 


9.4 


3.4 


0.37 


0.00 


Pchl3084. 


.3439 




Lipase 


16.0 


4.8 


0.3 


0.01 


Pchl3084. 


.3943 




Lipase 


25.9 


11.4 


0.44 


0.01 


Pchl3084. 


.0883 




Phospholipase 


50.8 


10.6 


0.21 


0.02 


Pyrrol nitrin 
















Pchl3084. 


.3143 


prnD 


Aminopyrrolnitrin oxidase 


6.94 


0.57 


0.08 


0.04 


Pchl3084. 


.3144 


prnC 


Halogenase 


12.53 


0.85 


0.07 


0.02 


Pchl3084. 


.3145 


prnB 


Pyrrolnitrin biosynthesis 


6.88 


1.23 


0.18 


0.00 


Pchl3084. 


.3146 


prnA 


Tryptophan halogenase 


16.45 


1.83 


0.11 


0.01 


Chitinase 
















Pchl3084. 


.2020 




Chitin-binding domain protein 


36.06 


4.28 


0.12 


0.03 


Pchl3084. 


.2021 




Chitinase 


31.17 


3.46 


0.11 


0.01 


Pchl3084. 


.3180 




Chitinase 


57.95 


1.83 


0.03 


0.02 


Pchl3084. 


.3181 




Chitin-binding domain protein 


51.78 


3.44 


0.07 


0.01 


Iron uptake 
















Pchl3084_3935 


pvdA 


L-ornithine 5-monooxygenase 


2.32 


5.11 


2.20 


0.01 


Pchl3084. 


.3936 


fpvl 


Sigma-70 factor 


12.77 


23.78 


1.86 


0.00 


Pchl3084. 


.0320 




TonB-dependent outer membrane receptor 


1.71 


6.19 


3.61 


0.02 


Pchl3084. 


.0855 




TonB-dependent outer membrane receptor 


23.45 


85.19 


3.63 


0.00 


Pchl3084. 


.2518 




TonB-dependent outer membrane receptor 


2.71 


7.17 


2.65 


0.00 


Pchl3084. 


.3174 




TonB-dependent outer membrane receptor 


8.26 


16.77 


2.03 


0.02 


Pchl3084. 


.3659 




TonB-dependent outer membrane receptor 


3.05 


6.43 


2.11 


0.00 


Pchl3084. 


.3848 




TonB-dependent outer membrane receptor 


7.11 


15.08 


2.12 


0.00 


Pchl3084. 


.5801 


ExbB 


TonB system transport protein 


55.52 


116.65 


2.10 


0.04 



Differential transcript abundance of 23 genes annotated 
as being involved in signal transduction or transcription 
regulation suggest that these genes also were regulated by 
GacA (Table S6). Among them, transcripts of 21 genes 
were underrepresented in the gacA mutant compared with 
WT. Most of these genes are functionally uncharacterized 
suggesting there remains a large gap in our understanding 
of P. chlororaphis gene regulation by the Gac system. 

Transcripts of genes involved in iron acquisition 
are in higher abundance in the gac mutant 
relative to WT: potential contribution to the 
Gac phenotype 

Hyper-fluorescence is a consistent phenotype of gac 
mutants among fluorescent pseudomonads, including 
30-84 (Chancey et al. 2002). It was shown previously that 
the Gac system negatively controls a large number of 
genes involved in iron acquisition (Hassan et al, 2010). 
Transcripts of nine genes annotated to be involved in iron 



acquisition were two- to threefold more abundant in the 
gacA mutant than the WT (Table 2). Six of them were 
annotated as TonB-dependent outer membrane receptors 
potentially involved in the uptake of ferric pyoverdine 
siderophores. According to current siderophore uptake 
models, when ferric pyoverdine is recognized by its cog- 
nate TonB outer membrane receptor, it is transported 
into the periplasmic space via a two-gate mechanism 
(Schalk 2008). TonB and proton motive force provide the 
energy for siderophore uptake and cause the release of 
the iron, which is subsequently transported into the cyto- 
plasm through an ABC transporter (Schalk 2008). Consis- 
tent with increased iron uptake, transcripts of pvdA 
(involved in fluorescent pyoverdine biosynthesis), exbB (a 
cytoplasmic membrane protein), and fpvl (an RNA poly- 
merase sigma-70 factor that recruits RNA polymerase for 
expression of iron transport genes) were overrepresented 
in the gacA mutant compared with WT. Previously it was 
shown that insertion of Tn5 into the pvdA locus of 30-84 
resulted in the loss of fluorescence on low iron media 
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(data not shown), linking fluorescence with pyoverdine 
synthesis. These results are consistent with overproduc- 
tion of pyoverdine being partially responsible for the 
hyper-fluorescent phenotype typical of gac mutants of 
Pseudomonas and suggest a possible benefit to having 
Gac~ phenotypes in mixed populations. 

Transcripts of genes involved in preparation for 
exponential growth are in higher abundance in 
the gac mutant relative to WT: GacA regulation of 
sigma factors and genes involved in translation 

Protein synthesis is a highly coordinated process that 
requires ribosomes composed of multiple proteins and 
ribosomal RNAs, mRNA, tRNAs, initiation and elonga- 
tion factors, and many other genes. One of the most 
dramatic transcriptional consequences of loss of gacA in 
P. chlororaphis 30-84 was the overexpression (two- to 
fourfold) relative to the WT of 50 genes involved in 
translation, including those encoding most of the large 
and small subunit ribosomal proteins, the translation ini- 
tiation factor IF-3, and the three translation elongation 
factors (Table S7). Because the ribosomal protein tran- 
scripts were increased in the gacA mutant, it is logical to 
speculate that rRNAs also were increased. However, due 
to the rRNA depletion steps, the RNA-seq only detected 
small amounts of rRNA transcripts. Therefore, we per- 
formed qPCR to quantify the relative abundance of the 
16S rRNA in the gacA mutant compared with the WT. 
Relative to WT, the abundance of 16S rRNA was signifi- 
cantly higher (1.8 ± 0.5-fold) in the gacA mutant, sug- 
gesting higher expression of ribosomal RNA. In addition, 
tRNA genes for transferring aspartate, glutamine, lysine, 
isoleucine, and alanine were overexpressed in the gacA 
mutant compared with WT; whereas threonine and gly- 
cine tRNAs were underexpressed (Table S7). A unique 
consequence of these findings was that due to the consis- 
tency of rpoD expression relative to 16S rRNA abun- 
dance, rpoD was used as the reference for all qPCR 
analyses comparing gene expression in the gacA mutant 
and WT. 

In previous study with 30-84 and other pseudomonads, 
it was found that mutations in gacA or gacS resulted in 
earlier transition of cells from lag to exponential phase 
growth, and this was suggested as an explanation for the 
competitive advantage of gac mutants in liquid broth 
cultures (Duffy and Defago 2000; Chancey et al. 2002; 
van den Broek et al. 2005; Driscoll et al., 2011). Although 
the precise mechanism(s) that facilitate gac mutants 
entering exponential growth phase earlier than the WT is 
unclear, we speculate that the reduction in rpoS expres- 
sion in gac mutants relative to WT may contribute to this 
effect. Differential regulation of the stationary phase 



sigma factor RpoS, but not the housekeeping sigma factor 
RpoD in a gacA mutant is particularly interesting given 
that RpoS competes with RpoD for binding to RNA poly- 
merase core enzymes (Farewell et al. 1998; Battesti et al. 
2011). This competition, in part, determines overall pat- 
terns of gene expression in the cell. Most research on 
RpoS has focused on its role in the transition from expo- 
nential to stationary phase growth. However, recently 
Rolfe et al. (2012) showed that during the transition from 
lag to exponential phase growth rpoS expression increased 
transiently, presumably to enable the cell to respond to 
oxidative challenges due to the introduction of the 
stationary phase cells into oxygenated fresh medium. In 
contrast, levels of RpoD increased twofold during this 
transition period. A reduction in levels of RpoS in gacA 
mutants would result in higher levels of expression of 
RpoD-regulated genes, which are involved in primary 
growth. Increased levels of transcription of translational 
machinery, including ribosomal proteins, rRNAs, tRNAs, 
and accessory proteins in the gac mutant as reported here 
also would favor earlier entry into exponential growth. 

The role of GacA in bacterial motility 

The Gac system has been shown to differentially control 
motility in various bacteria. For example, it positively 
contributes to motility in E. coli K-12 (Wei et al. 2001), 
P. protegens Pf-5 (Hassan et al. 2010), and P. syringae pv. 
tabaci (Marutani et al. 2008), but negatively in P. fluores- 
ceins Fl 13 (Martinez-Granero et al. 2012) and Erwinia 
amylovora Eall89 (Zhao et al. 2009). In both Pf-5 and 
P. syringae pv tabaci, the Gac system positively influences 
the level of transcription of the flagellar genes (Marutani 
et al. 2008; Hassan et al. 2010). Consistent with the 
potential role of GacA in enhanced bacterial motility, col- 
onies of 30-84 gac mutants typically appear larger on 
standard culture plates. Bacterial motility was assessed for 
strain WT, AgacS, AgacA, and 30-84ZN by inoculating 
bacterial cells onto motility plates (0.25% agar) and mea- 
suring the diameter of the area colonized by bacterial cells 
at 24 and 48 h as described previously (Wang et al. 
2011). The Gac mutants exhibited enhanced motility 
compared with that of the WT strain 12 and 24 h follow- 
ing inoculation (Fig. S1A and B). In contrast, phenazine 
deficient derivative 30-84ZN (phzBr.lacZ) exhibited a sim- 
ilar motility pattern to the WT suggesting that the 
increased motility in the gac mutant was not due to lack 
of phenazine production. Diameters of colonies for the 
WT, AgacS, AgacA, and 30-84ZN at 48 h were about 
4.8 ± 0.1, 6.4 ± 0.3, 6.5 ± 0.4, and 4.8 ± 0.3 cm (lon- 
gest dimension from the inoculation point), respectively. 
These results indicate that the Gac system negatively regu- 
lates bacterial motility in 30-84. 



2013 The Authors. MicrobiologyOpen published by John Wiley & Sons Ltd 



519 



The Gac System in P. chlororaphis 

It was demonstrated previously that swimming motility 
in P. fluorescens F113 is under negative control by Gac 
and that this downregulation occurs through the repres- 
sion of the fiagella master regulatory gene fleQ (Martinez- 
Granero et al. 2012). One of the factors affecting motility 
is fiagella production. Moreover, a P. chlororaphis 30-84 
fliM mutant is nonmotile (L. S. Pierson et al., unpubl. 
data). However, RNA-seq analyses revealed that the tran- 
script abundance of flagellar genes such as fliA, fliM, and 
fliC were not significantly altered in the gacA mutant 
(Fig. SIC), and these data were verified by qPCR (Fig. 
SID). The results suggest that the impact of GacA on 
motility may be secondary to fiagella structure, possibly 
related to flagellar rotation. Another possibility is that 
increased swimming motility is related to enhanced 
surfactant production, as observed in other bacterial spe- 
cies (Wilf et al., 2012; Daniels et al. 2004). However, 
known surfactant biosynthetic genes such as the rhlA or 
homologs were not detected in P. chlororaphis 30-84 
genome. 

Conclusions 

The role of the ncRNAs rsmY and rsmZ in the Pseudomo- 
nas Gac/Rsm signal transduction pathway was first 
detected by their binding capacity to the regulatory pro- 
tein RsmA and subsequently RsmE in P. fluorescens 
CHA0, where their interaction contributed to the regula- 
tion of HCN, exoprotease, and 2,4-diacetylphloroglucinol 
biosynthesis important for the biological control capabil- 
ity of the strain (Heeb et al. 2002; Reimmann et al. 
2005). In our study, the functionality of the components 
of the Gac/Rsm signal transduction pathway in the regu- 
lation of phenazines and other genes important for 
biological control in strain 30-84 was examined. Only 
rsmZ and rsmE are linked to phenazine production, how- 
ever, further RNA-protein binding experiments are 
needed to verify whether direct interactions between 
RsmE and the ncRNA rsmZ are involved in the regulation 
of phenazine biosynthesis. Interestingly, in P. fluorescens 
CHA0 rsmE expression was shown to be regulated nega- 
tively by RsmA (Reimmann et al. 2005). What role RsmA 
plays in Gac-mediated regulation in 30-84 has yet to be 
determined. In P. fluorescens CHA0 and in P. aeruginosa, 
the regulatory effect of gacA was attributed exclusively to 
the direct control of only two ncRNAs rsmY and rsmZ 
(Heeb et al. 2002; Brencic and Lory 2009). In P. aerugin- 
osa, single rsmY and rsmZ mutants synthesize intermedi- 
ate levels of phenazines, whereas the double rsmYlrsmZ 
mutant produces less (Kay et al. 2006). Although constit- 
utive expression of rsmZ completely restored phenazine, 
HCN, and exoprotease production in a 30-84 gacA 
mutant, we cannot exclude the possibility that rsmY also 
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partially contributes to the production of these key bio- 
logical control compounds. Future study will be directed 
toward identifying the potential roles of the three ncRNAs 
in strain 30-84. 

In most Pseudomonas species, phenazine production is 
controlled by combinations of conserved regulatory sys- 
tems integrated into sensory networks potentially respon- 
sive to environmental, nutritional, population, and 
metabolic inputs (reviewed in Pierson and Pierson 2010). 
Transcriptomic analysis revealed that the expression of 
regulatory genes known to be involved in phenazine 
production (including phzR and phzl, pip, and rpoS) as 
well as other genes annotated as being involved in 
trascriptional regulation were affected by a gacA in muta- 
tion 30-84. Our genetic analysis demonstrated that 
whereas a gacA mutation could be bypassed by either 
mutation of rsmE or constitutive expression of rsmZ, it 
could not be rescued by expression of quorum sensing 
(e.g., constitutive expression of phzR with the addition of 
AHLs) or other regulatory genes. These data provide 
evidence that the phenazine biosynthesis genes are direct 
targets of Gac/Rsm regulation and highlight the strength 
of gene deletion/overexpression analysis in combination 
with transcriptomic profiling as an approach to study 
complex regulatory mechanisms. These results also illus- 
trate the multiple-level regulatory role of the Gac/Rsm 
system in the control of phenazines and suggest that the 
far-reaching influences of GacA on the transcriptome are 
likely mediated through intermediate transcriptional 
regulatory genes. 

Our transcriptional profiling revealed that GacA regu- 
lates a wide array of biological functions in P. chlororaphis 
30-84. In this study, transcriptomic analysis identified 771 
genes that were differentially expressed in the GacA 
mutant compared with WT 30-84, with 551 genes down- 
regulated in the gacA mutant under our experimental 
conditions. Consistent with previous observations that gac 
mutants are phenazine, HCN, exoprotease, and lipase 
deficient, the transcript abundances of genes that code for 
these phenotypes were reduced in the gacA mutant 
compared with WT 30-84. Additionally, the transcript 
abundances of 30-84 genes annotated to encode pyrrolni- 
trin biosynthesis and chitin degradation also were 
reduced. These results are consistent with previous find- 
ings demonstrating that GacS/GacA functions as a posi- 
tive regulator of secondary metabolites that contribute to 
biological control activity in Pseudomonas. Other genes 
positively regulated by GacA in strain 30-84 that could be 
of potential importance in rhizosphere survival, plant- 
microbe interaction, and biological control include genes 
annotated as coding for components of a T6SS. Beyond 
genes involved in well-defined functional processes, the 
gacA regulon included a gene cluster that is involved in 
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the utilization of polyhydroxyalkanoates, important 
storage molecules in many pseudomonads. Further stud- 
ies will be required to determine whether either strain 
contains a functional T6SS and whether strain 30-84 is 
capable of utilizing polyhydroxyalkanoates. Future studies 
also will be directed at how Gac-mediated regulation of 
these potentially beneficial genes is altered under rhizo- 
sphere conditions. 

Many Pseudomonas species identified for biological 
control demonstrate phenotypic variation resulting from 
spontaneous mutation in gacS or gacA. However, pheno- 
typic variation derived from conserved mutation of spe- 
cific regulatory genes is not limited to Pseudomonas, and 
has been observed in many Gram - and Gram + species 
resulting from mutations in gac or other regulatory genes 
(van den Broek et al. 2005). The prevalence of Pseudomo- 
nas gac mutants in laboratory and natural environments 
suggests that there is some benefit to the gac mutant 
phenotype and thus selection at some level to maintain 
phenotypic variation in the species. Previous study with 
strain 30-84 provides evidence that the presence of the 
gac mutants in WT-Gac mixed populations in biofilms 
and on plant roots can be beneficial (Chancey et al. 2002; 
Driscoll et al, 2011). Thus, of particular interest in this 
study was the identification of genes that are expressed at 
a higher level in 30-84 gacA mutants. Our transcriptomic 
analysis showed that 220 genes that were overexpressed in 
the gacA mutant under our experimental conditions. The 
most obvious group, genes involved in iron uptake, may 
increase iron availability especially in the biofilm context. 
Increased protein translation is another possible benefit, 
as the shorter transition to exponential growth for gac 
mutants allows them to quickly establish a population. It 
is worth noting that 52 upregulated genes in the gacA 
mutant are hypothetical or unknown genes, including the 
highest upregulated gene (Pchl30-84_3493, 26.2-fold 
increase) (Table S3). Functional characterizations of these 
highly expressed genes are currently underway. 
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Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Figure SI. (A) Migration phenotypes of WT, gac mutants 
and 30-84ZN at 48 h after incubation. (B) Comparison of 
the movement distance of WT, gac mutants and 
30-84ZN. Cells were inoculated on the surface of a motil- 
ity plate containing 0.25% agar as described previously 
(Wang 2011). Plates were incubated at 28°C for 48 h, 
during which the displacement diameter to the outermost 
edge of the movement area was measured. The experi- 
ment was performed at least three times in triplicate. (C) 
Expression levels of fliA, fliC, and fliM in WT and a gacA 
mutant. RPKM values representing levels of gene expres- 
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sion were calculated from RNA-seq analysis and log 
transformed. Data points represent means of three repli- 
cates ± standard deviations. (D) Verification of RNA-seq 
data with qPCR. Bacterial strains were grown overnight 
in LB broth and reinoculated in 5 mL AB medium + 2% 
casamino acid. Relative expression of fliA, fliC, and fliM 
genes, normalized to the expression of rpoD gene, was 
determined by qPCR after ~18 h growth (OD 600 at 1.2). 
Table SI. Bacterial strains and plasmids used in this 
study. 

Table S2. Oligonucleotides used for gene cloning and 
qPCR. 

Table S3. Differentially expressed genes in the gacA 
mutant compared with the WT strain. 
Table S4. Mean transcript abundance and ratio of abun- 
dances (AgacA/WT) of type VI secretion system (T6SS) 
genes in the gacA mutant compared with the WT. 
Table S5. Mean transcript abundance and ratio of abun- 
dances (AgacA/WT) of polyhydroxyalkanoate (PHA) bio- 
synthetic genes in the gacA mutant compared with the 
WT. 

Table S6. Mean transcript abundance and ratio of abun- 
dances (AgacA/WT) of regulatory genes in the gacA 
mutant compared with the WT. 

Table S7. Mean transcript abundance and ratio of abun- 
dances (AgacA/WT) of genes involved in protein metabo- 
lism in the gacA mutant compared with the WT. 
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